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Abstract. We present the full characterization of phase-randomized or phase- 
averaged coherent states, a class of states exploited in communication channels and 
in decoy state-based quantum key distribution protocols. In particular, we report 
on the suitable formalism to analytically describe the main features of this class of 
states and on their experimental investigation, that results in agreement with theory. 
We also show the results we obtained by manipulating the phase-averaged coherent 
states with linear optical elements and testify their good quality by employing some 
non-Gaussianity measures and the concept of mutual information. 
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1. Introduction 

Laser radiation, which can be described in terms of coherent states, plays a relevant role 
in practical communication schemes. One of the main advantages of coherent states 
over more exotic quantum states, such as squeezed ones, is that they can propagate 
over long distances, only suffering attenuation and without altering their fundamental 
properties. A coherent state is characterized by a Poissonian photon-number statistics 
and a well-defined optical phase. Thus, one can easily implement phase-shifted keyed 
communication in which the logical information (the bit) is encoded in two coherent 
states with the same amplitude and a vr-difference in phase. Nevertheless, this kind 
of communication channel lacks security. Remarkably, very recently quantum key 
distribution involving coherent states and decoy states has been realized and it has 
been pointed out that phase-averaged coherent states may enhance the security of the 
channel [H |2l |3] . In this case, the high degree of accuracy in the phase randomization 
process is one of the main requirements. 

By contrast to coherent states, which are described by Gaussian Wigner functions, 
phase-averaged coherent states clearly exhibit non-Gaussian features [4]. Thus, the 
systematic study of the nature of these states and the possibility to manipulate them 
can be considered of real interest in enhancing the performances of the communication 
protocols in which they are employed [5]. 

In this paper we investigate the main features of phase-averaged coherent states and 
report on their fully experimental characterization by addressing the measurement of the 
photon-number statistics and the reconstruction of the Wigner function. Furthermore, 
we perform basic manipulation experiments by means of linear optical elements, in order 
to assess the usefulness of these states for communication and information processing. 
The detection of such states is performed in the mesoscopic photon-number regime by 
means of hybrid photodetectors. 

The plan of the paper is as follows. In section |2] we recall several properties of phase- 
averaged coherent states like photon-number statistics and purity. We also present the 
experimental scheme used for the generation, characterization and manipulation of such 
states. Section [3] is devoted to their basic description using the characteristic functions 
and the corresponding quasiprobability densities. In this section we also report on the 
strategy and realization of the experimental reconstruction of the Wigner function. Non- 
Gaussianity of the phase-averaged coherent states and its experimental measurement 
are addressed in section |H Section [5] investigates the linear operations with phase- 
averaged coherent states performed by a beam splitter. Here we give a complete analytic 
description of the output one-mode reduced states. These superpositions turn out to 
be Fock-diagonal and are interesting for quantum information processing. We find a 
good agreement between theory and our experimental results for what concerns non- 
Gaussianity and mutual information of the beam-splitter output states. Our concluding 
remarks are drawn in section [6l 
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2. Quantum description of phase-averaged coherent states 



A single-mode phase-randomized or phase-averaged coherent state (PHAV) is obtained 
by randomizing the phase of a coherent state: 



exp { --I 



n=0 



(1) 



with (3 



where 



j\^^-y pfjy^Y g ig diagonal in the photon-number basis, namely: 
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is a Poissonian distribution. Therefore, randomizing the phase of a coherent state does 
not change its photon-number distribution [see ([3])]. Moreover, due to the diagonal 
structure of its density matrix, the statistical properties of a PHAV can be fully described 
by the photon- number distribution. Indeed, while is a pure state, the degree of purity 
hIq] of the PHAV q can be directly evaluated by the photon- number distribution and is 
given by: 



M^] = E^L = exp(-2|/3p)/o(2|/3p), 



(4) 



n=0 



where Iq{z) is the modified zeroth-order Bessel function of the first kind. The purity ^[q] 
is a strictly decreasing function of the mean number of photons {p)a) = Tr[^a''"a] = 
(a is the annihilation operator and [a, a''] = I). 

From the experimental point of view, we obtained this class of states by sending the 
second-harmonics pulses of a mode-locked Nd:YLF laser amplified at 500 Hz (High-Q 
Laser Production) to a mirror mounted on a piezo-electric movement (see figured]). The 
displacement of the piezoelectric movement, which is controlled by a function generator, 
is operated at a frequency of 100 Hz and covers a 12 fim span [6]. 




Figure 1. Experimental setup. F^: variable neutral density filter; BS: 50:50 beam 
splitter; Pz: piezoelectric movement; MF: multimode fiber (600 fim core); HPD: hybrid 
photodetector. 



In figure |2] we show the detected photon distributions of three PHAVs at different 
energy values, obtained by using a direct detection scheme employing a hybrid 
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photodetector (HPD, R10467U-40, maximum quantum efficiency ~ 0.5 at 500 nm, 
Hamamatsu) characterized by a partial photon-counting capabihty and a hnear response 
up to 100 photons [HIH]- In the same figure, we also show the corresponding theoretical 
photon-number statistics for detected photons [the photocount distribution is simply 
obtained by using ([3]) and replacing by M = r]\(3\'^, where rj is the overall 
quantum efficiency]. The good agreement between experimental data and theory can 




Figure 2. Detected-photon distribution of PHAV state for three different mean 
values. Coloured dots: experimental data, lines: theoretical expectations. The purity 
is n[g] = 0.13 (red plot), iJ.[g] = 0.20 (blue plot) and n[g] = 0.39 (black plot). 

be quantified by calculating the fidelity (see F values reported in figure [2]): F = 
X]m=o \/Pth{TTT')P{^), in which Pth{m) and P{m) are the theoretical and experimental 
distributions, respectively, and the sum is extended up to the maximum detected-photon 
number fh above which both P^-h{m) and P{m) become negligible. 

3. Basic characterization 

We insert ([2]) into the well-known definition of the characteristic functions (CFs): 



to write: 



X(A; s) = exp ( ||Apj Tr[^^(A)], (-1 < s < 1) (5) 

/ 1 „ s OO 

x(A;.)=exp — IAH ^nn lJ\X\') . (6) 

^ ^ n=0 

The series (E]) is obtained by substitution of the diagonal matrix elements of the 
displacement operator D{X) = exp(Aa^ — A*d) in the photon-number basis, namely: 

HD(A)|n) = exp(-i|Ap) {\X\') , (7) 
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where Ln{z) is a Laguerre polynomial. By using the generating function of the Laguerre 
polynomials [9l |T0] : 



^L„(z)^ = eVo(2 



n=0 



where Jq{z) is a Bessel function of the first kind, we get: 

1 - s, 



X(A; s) = exp 



An Jo(2|/3||A|). 



(8) 



(9) 



The Fourier transform of the CFs in (Q gives us the set of quasiprobability densities 

mm- 



W{z; s) = - f d^A exp (^A* - z*X) x(A; s), 
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which exphcitly read P [TO] : 
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The P function retrieved for s = 1 by employing asymptotic expansions of Bessel 
functions is expressed in terms of Dirac's 6 distribution: 

P{z) = W{z;l) 



Is 

71 



27r| 



(13) 



and the last equality follows from the properties of the 6 distribution. Equation (fT3l) 
shows us that the mixed state obtained by averaging over the phase of a pure coherent 
state preserves the important feature of being at the classicality threshold (remember 
that the coherent states are the only pure states at the classicality threshold). 
Finally, for s = we get the Wigner function: 



W{z) = 2 exp [-2 ( 



+ z 



(14) 



which is positive everywhere in the phase space. Recently, states with positive Wigner 
functions have become interesting for efficient classical simulation of a broad class of 
quantum optics experiments. In [T3l [H] a protocol for classical simulations using non- 
Gaussian states with positive Wigner function was presented (see also the more recent 
[15]). Note that the Wigner function in ( HM is not Gaussian, a feature that becomes 
evident from the plot of the function that shows a dip at the origin of the phase space 
(see figure [3]). 

In order to experimentally reconstruct the Wigner function of PHAVs we adopted 
the same strategy presented in [T6j and based on the measurements of the statistics of 
the state under investigation mixed at a beam splitter (BS) with a coherent probe field 
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whose amplitude and phase can be continuously changed [171 dSl [ISj. As the PHAV 
is a diagonal state, its Wigner function is phase insensitive, i.e. it exhibits a rotational 
invariance about the origin of the phase space. For this reason, in figure [3] we show 




Figure 3. Experimental reconstruction of a section of the Wigner function of a PHAV 
with |/3p — 1.97 and ^ — 0.999. Blue dots: experimental data, orange mesh: theoretical 
expectation. 



the experimental data (blue dots) corresponding to a section of the Wigner function 
superimposed to the theoretical surface (orange mesh): 



W \ ia ]=W[ \ia exp 



1-e 



(15) 



^ being the overall (spatial and temporal) overlap between the probe field and PHAV 
[T6] . In ( !T5|) . is now the mean number of photons we measured, which includes the 
quantum efficiency. In fact, it is worth noting that for classical states the functional form 
of the Wigner function is preserved also in the presence of losses and its expression, given 
in terms of detected photons, reads W{a) = ^Y^m=oi~^)"^Prna^ where represent 
the detected-photon-number distributions of the state to be measured displaced by the 
probe field [16]. 



4. Non-Gaussianity of PHAVs 

In order to quantify the non-Gaussanity of PHAVs, here we compare three different 
measures of non-Gaussianity recently introduced [201 1211 [221 [23] and completely 
determined by the density matrix in ([2]) and ([3]). All the three measures compare 
the properties of the state under investigation g with that of a Gaussian reference state, 
a, having the same mean value and covariance matrix as g. In the case of PHAVs, the 
reference Gaussian state is a thermal state with mean occupancy |/3p. 
The first measure is based on the Hilbert-Schmidt distance: 
^ Dls[0^ ^] - 2K[g,a] 

H[g\ 2/i[^J 
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where id is the purity of the state and K[g, a] =Tr[^(j] [20]. The Hilbert-Schmidt distance 
can be analytically calculated by using the purity (jl]), the degree of purity of the reference 
thermal state = (2|/3p + 1)"^ and the expression of n[g,a]: 

1 f \ 



K[g,a] 



exp 



(17) 



The second measure is the relative entropy of non-Gaussianity defined as: 

eB[g] = Sia)-Sig), (18) 

where S{g) = — Tr[^ In g] is the von Neumann entropy of the state g [21]. For all diagonal 
states, we have S{g) = —J2n Qnn^^ gnn, where in the present case gnn is given in (j3]), 
and Sia) = (|/3|2 + 1) In (|/3p + 1) - l/^pln 

The last measure we study has been recently introduced in [23] and is based on the 
quantum fidelity, namely: 

ec[g] = l-'^J'ig,a), (19) 

where: 



J^{g,a) 



Tr 



' g(y \J g 



(20) 



is the Uhlmann fidelity [21]. This measure can readily be evaluated for Fock-diagonal 
states since they commute with their reference thermal states. 

It is worth noting that the evaluation of the considered measures only requires 
quantities that can be experimentally accessed by direct detection, as they can be 
expressed in terms of photon- number distributions [26] . 



««3 




10 12 



Figure 4. Comparison between ttie tliree measures of non-Gaussianity in the case 
of PHAVs as functions of the mean number of photons. Open squared symbols: 
experimental data, dots: theoretical expectations. 

Figure m shows the behaviour of the three non-Gaussianity measures as functions of 
the average number M of detected photons in the case of PHAVs. It is evident that the 
behaviours of the three measures are very similar to each other except for the absolute 
values [27] . 
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Figure 5. Mutual information between the two outputs of the BS at which a PHAV 
with energy Mt is divided. Open squared symbols: experimental data, red dots: 
theory. The error bars are smaller than the symbol size. 



5. Advanced characterization and manipulation 

When a coherent state is mixed with the vacuum at a BS with transmissivity r, the 
two emerging beams are excited in the product state |\/r/3) ® |a/1 — r/3) and thus are 
uncorrelated. Nevertheless, when we consider a PHAV as the input state, a correlation 
arises at the two outputs, even if intensity correlations still vanish [28]. The total amount 
of correlation of the output bipartite state Qu can be evaluated in terms of the mutual 
information (MI): 



where Qk = TihiQu] {h, k = 1,2 and k ^ h) are the output reduced states and S{q) is the 
von Neumann entropy. In figure [5] we plot the experimental data (open squared symbols) 
and the theoretical predictions (red dots) of the MI as a function of the energy of the 
input PHAV. We experimentally measured this parameter by using a scheme involving 
two hybrid photodetectors to simultaneously detect the light at the two outputs of the 
BS [29], as shown in figure [H Since from the experimental point of view it is not possible 
to measure both the input and output states simultaneously, in order to assess the last 
term in fl2ip we assumed that the input state was a PHAV with energy equal to the 
sum of the two output channels (losses at the BS are negligible). We also notice that in 
figure the slight discrepancy between experiment and theory appearing at increasing 
values of the input energy can be due to some saturation effect of the acquisition chain. 

We have also investigated another interesting state obtained by the interference 
of two PHAVs (see figure [1]): we will refer to this state, which is still diagonal in 
the photon-number basis, as 2-PHAV [26]. The 2-PHAV can find useful applications 
in passive decoy state quantum key distribution [5|. To describe the 2-PHAV state, 
we start observing that when two uncorrelated field modes described by a product 
CF xo(Ai,A2;s) = Xoi(Ai; s) ^02(^2; s) are mixed at a BS with transmissivity r, the 



(21) 
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output two-mode CF may be written as follows x(Ai, A2; s) = Xoi(Ci) s) Xo2{C2', s) with 
(1 = Aiv/^ - AsTT^, C2 = Aav^ + XiVT^ [1]. Therefore, the CF of the 2-PHAV 
state, obtained by taking only one of the output modes, can be formally written as (the 
CF of the other mode can be obtained by replacing r with 1 — r): 

X(A; s) = xoi (Av^; s) xm [\VT^; s) , (22) 

which follows from the partial-trace rule in the reciprocal phase space. f l22|) gives the 
following multiplication rule for the input states of the type (j9]): 

1 - s, 



X(A; s) = exp ^l^l j^o [2\Pi\\\\V^) -h [2\h\\\\VT^) , (23) 

Pi and are the coherent amplitudes of the corresponding interfering PHAVs. Note 
that x(A; s) = x(|A|; s), as expected for phase-insensitive states. 

In order to obtain the photon statistics of the 2-PHAV, we can follow two strategies. 
On the one hand, we can use the expansion of the density operator in terms of 
displacement operators [TT] : 

g=- [ d2Ax(A;0)D(-A). (24) 



TT JC 

Therefore, the matrix elements of a phase-insensitive single-mode state described by the 
CF x(|A|) = x(A; 0) may be written as: 



Qlr, 



2 Si^ d|A| |A| exp (-^|Ap) x(|A|) L„ (|Ap) , (25) 

where we used ((Tj) and performed the integration over the polar angle, thus being left 
with an integral over | A | . 

On the other hand, we can exploit high-order correlation functions via the series 

mm-- 

1 00 / 1 \ k—n 

.... = ;^i:^((«*)'a'>. (26) 

Using a generating-function method, we are able to derive the /cth-order correlation 
functions in terms of Legendre polynomials Pk{x): 

((at)'=a'=) = + |/32p(l - r)f P, Q , (27) 

where u is the ratio: 

u 





V - 1/32 


^(l-r)l 






2(1 -r) 



The fcth-order normalized correlation functions are then: 

,<..(0) = i(^ = „'=pJiV (29) 



{a}a)^ \u, 

Now, thanks to the relationship ( l26l) we obtain the density matrix elements as series 
expansions involving Legendre polynomials: 
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In particular, if r = 1/2 (balanced BS) and = |/32| = (identical inputs), we find 
a simple result for the degree of coherence fl29l) . namely: 

9^'\0) = ^^^j^ > I, ik>l), (31) 

which indicates a super-Poissonian photon statistics. The photon-number distribution 
is obtained after some algebra via (126|) as: 



: ^^1^ exp (-2|/3p) (1; n + 1; 21/3^) , (32) 
where iFi(g;p;2) is a confluent hypergeometric function pi [30]. 



[n\ 




Figure 6. "Detected-photon distribution of a balanced 2-PHAV state for three 
different mean values. Coloured dots: experimental data, lines: theoretical 
expectations (|32|) . The purity is ii[g] = 0.09 (red plot), ij[q] = 0.13 (blue plot) and 
fi[g\ = 0.21 (black plot). 

In figure |6] we have plotted the photon-number distribution (32) for different energy 
values (coloured lines). We remark the good agreement between experimental data 
(coloured dots) and theory predictions confirmed also by the high values of the fidelity. 

Finally, we have obtained the set of quasiprobability densities associated with a 
2-PIIAV as having the following expansion: 



W{a; s) = 2 exp 



2\a\ 



where Pk{z) are Legendre polynomials and u is defined in (1281) . In particular, the 
quasiprobability densities of the balanced state (|32|) have a simpler form due to the 
explicit correlation functions ( !3TI) . We get: 

.... X / 2|«n - (2A;-1)!! / 2|/3P\\ f2\a\^\ 
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W(z) 



Figure 7. Left: Experimental reconstruction of a section of the Wigner function of a 
PHAV with = 1.97 and S, = 0.999. Blue dots: experimental data, orange mesh: 
theoretical expectation. Right: Experimental reconstruction of a section of the Wigner 
function of a balanced 2-PHAV with = 1.03, |/32p = 0.91, = 0.95 and = 1- 
Red dots: experimental data, blue mesh: theoretical expectation. 



In figure [7]we report a section of tlie pliase- insensitive Wigner function of a 2-PHAV 
(on the right) obtained by the interference at a balanced BS of two identical PHAVs 
(whose section of Wigner function is shown on the left). Even in the case of 2-PHAV, 
the experimental data (red dots) are well superimposed to the theoretical surface (blue 
mesh) : 

W2^PRAY{^/^a) = W^2-PHAv(\/ep«) exp [ - \a\^l - - (|/3i| + - (35) 

where describes the overall overlap between the probe and the 2-PHAV and .^s the 
overall overlap between the two components of the 2-PHAV. It is evident that the single 
PHAV has a dip at the origin of the phase space, whereas the 2-PHAV has a peak. This 
difference results in a reduction of non-Gaussianity of the 2-PHAV with respect to that of 
a single PHAV at fixed energy, as testified by the non-Gaussianity measures introduced 
above [261 [27]. To stress this result, in figure [H] we show the behaviour of the three 
measures as functions of the energy values M for a balanced 2-PHAV: we can notice 
that the absolute values of Ek, k = A, B,C are smaller than the ones we obtained in the 
case of a single PHAV. Moreover, in figure |9] we plot the same measures as functions of 
the balancing between the two components of the 2-PHAV at fixed energy value M. As 
one may expect, the three measures monotonically decrease at increasing the balancing. 
In fact, the most unbalanced condition reduces to the case in which there is only a 
single PHAV, whereas the most balanced one corresponds to have a balanced 2-PHAV. 
As in the case of a single PHAV, when a 2-PHAV state is mixed with the vacuum at 
a 50:50 BS, the two outputs show a correlated nature, testified by the non-zero value 
of the mutual information. In figure [10] we report the MI between the two outputs of 
the BS at which a balanced 2-PHAV is divided as a function of the input energy value. 
Furthermore, figure [TT] shows the MI at fixed input energy of the 2-PHAV as a function 
of the ratio R = |/3i|/|/32| between the two single PHAVs used to generate the 2-PHAV 
state. 
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Figure 8. Comparison between tlie three measures of non-Gaussianity in the case 
of balanced 2-PHAVs as functions of the mean number of photons. Open squared 
symbols: experimental data, dots: theoretical expectations. 
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Figure 9. Comparison between the three measures of non-Gaussianity in the case of 
2-PHAVs at fixed total energy as functions of the balancing between the two single 
PHAVs. Open squared symbols: experimental data, dots: theoretical expectations. 



In both figures the accordance between the experimental data (open squared 
symbols), whose error bars are smaller than the symbol size, and the theoretical 
predictions (red dots) is good. As in the case of a single PHAV, in order to assess 
the last term in ( 12T|) we again assumed that the input state was a 2-PHAV with energy 
equal to the sum of the two output channels. 

6. Concluding remarks 

In conclusion, we have studied the main properties of PHAVs: we have developed 
an analytic description and verified the theoretical predictions by means of a direct 
detection scheme involving hybrid photodetectors. In detail, we have investigated the 
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Figure 10. Mutual information between the two outputs of the BS at which a balanced 
2-PHAV is divided as a function of the total input energy AIt- Open squared symbols: 
experimental data, red dots: theory. The error bars are smaller than the symbol size. 
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Figure 11. Mutual information between the two outputs of the BS at which a 2-PHAV 
at fixed input energy Mt is divided as a function of the ratio R = |/3i|/|/32| between 
the energies of the two PHAVs generating the 2-PHAV state. Open squared symbols: 
experimental data, red dots: theory. The error bars are smaller than the symbol size. 



detected photon-number distribution and the Wigner function that is non-Gaussian. 
Moreover, we have used three different non-Gaussianity measures, all based on quantities 
experimentally accessed by direct detection, to quantify the non-Gaussianity amount 
and have proven the consistency of the different approaches. Furthermore, we have 
manipulated PHAVs by means of linear optical elements and generated a new class of 
phase-randomized states, namely 2-PHAVs, obtained as superpositions of two PHAVs 
at a BS. The consistent experimental and theoretical results we have obtained in 
the characterization of both PHAVs and their superpositions 2-PHAVs reinforce the 
possibility of using them for applications to communication protocols. These classical 
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states appear to be robust, experimentally accessible and theoretically convenient. The 
investigation of the non-Gaussianity and correlations of some other output BS-states 
manipulated by conditional measurements is one of our present interests. 
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